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Radio metric data from the Pioneer 10/11, Galileo, and Ulysses spacecraft indicate an 
apparent anomalous, constant, acceleration acting on the spacecraft with a magnitude ~ 

8.5 x 10~ 8 cm/s 2 , directed towards the Sun. Two independent codes and physical strategies 
have been used to analyze the data. A number of potential causes have been ruled out. We 
discuss future kinematic tests and possible origins of the signal. 
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Exploration of the outer planets began with the launch 
of Pioneer 10 on 2 March 1972, JIJ . (Pioneer 11 followed 
on 5 April 1973.) After Jupiter and (for Pioneer 11) Sat¬ 
urn encounters, the two spacecraft followed hyperbolic 
orbits near the plane of the ecliptic to opposite sides of 
the solar system. Although Pioneer 10 is still transmit¬ 
ting, its mission officially ended on 31 March 1997 when 
it was at the distance of 67 Astronomical Units (AU) 
from the Sun. Pioneer ll’s radio system failed and co¬ 
herent Doppler signals were last received on 1 October 
1990, when the spacecraft was 30 AU away from the Sun. 

The Pioneer spacecraft are excellent for dynamical as¬ 
tronomy studies. Due to their spin-stabilization and their 
great distances, a minimum number of Earth-attitude re¬ 
orientation maneuvers are required. This permits precise 
acceleration estimations, to the level of 10~ 10 cm/s 2 (av¬ 
eraged over 5 days). Contrariwise, a Voyager-type space¬ 
craft is not well suited for a precise celestial mechanics 
experiment as its numerous attitude-control maneuvers 
overwhelm any small external acceleration. 

To obtain the S-band Doppler data from the Pioneer 
spacecraft, NASA used the Jet Propulsion Laboratory’s 
(JPL) Deep Space Network (DSN). This data was used 
in the two analyses described below to determine Pio¬ 
neer’s initial position, velocity and the magnitudes of the 
orientation maneuvers. The analyses were modelled to 
include the effects of planetary perturbations, radiation 
pressure, the interplanetary media, general relativity, and 
bias and drift in the range and Doppler. Planetary coor¬ 
dinates and the solar system masses were obtained using 
JPL’s Export Planetary Ephemeris DE200. Both analy¬ 
ses calculated Earth’s polar motion and its non-uniform 
rotation using the International Earth Rotation Service. 

Beginning in 1980, when at 20 AU the solar radiation 
pressure acceleration had decreased to<5xl0~ 8 cm/s 2 
[||, JPL’s Orbit Determination Program (ODP) analy¬ 
sis of unmodelled accelerations (at first with the faster- 


moving Pioneer 10) found that the biggest systematic 
error in the acceleration residuals is a constant bias of 
ap ~ (8 ± 3) x 10~ 8 cm/s 2 , directed toward the Sun [|j, 
to within the accuracy of the Pioneers’ antennae. As 
possible “perturbative forces” to explain this bias, we 
considered gravity from the Kuiper belt, gravity from 
the galaxy, spacecraft “gas leaks,” errors in the plane¬ 
tary ephemeris, and errors in the accepted values of the 
Earth’s orientation, precession, and nutation. None of 
these “forces” explained the apparent acceleration. Some 
were three orders of magnitude or more too small. 

Non-gravitational effects, such as solar radiation pres¬ 
sure and precessional attitude-control maneuvers, make 
small contributions to the apparent acceleration we have 
observed. The solar radiation pressure decreases as r~ 2 . 
As previously indicated for the Pioneers, at distances 
>10-15 AU it produces an acceleration that is much less 
than 8 x 10~ 8 cm/s 2 , directed away from the Sun. (The 
solar wind is roughly a factor of 100 smaller than this.) 

A possible systematic explanation of the residuals is 
non-isotropic thermal radiation. Pu 238 (half life of 87.74 
years) radioactive thermal generators (RTGs) power the 
Pioneers. At launch the RTGs delivered 160 W of elec¬ 
tric power. Power has decreased approximately linearly 
ever since. By 1997 a little less than 80 W were available. 
The excess power and the heat generated by the pluto¬ 
nium has been thermally radiated into space. The power 
needed for this to explain ap is ~ 85 W. There is almost 
that much available, but presumably the radiation was 
approximately isotropic. Further, if it were not, and was 
the cause of ap, this acceleration would have decreased 
with time. After 1980, no such (linearly decreasing) ac¬ 
celeration was observed. Another radiation source is the 
Pioneer radio beam. The power emitted from the an¬ 
tenna is 8 W. This implies a bias maximum of less than 
9% of ap, and in the opposite direction. (The influence 
of the bias is being investigated.) 
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We conclude, from the JPL-ODP analysis, that there 
is an unmodelled acceleration, ap, towards the Sun 
of (8.09 ± 0.20) x 10 -8 cm/s 2 for Pioneer 10 and of 
(8.56±0.15) x 10 -8 cm/s 2 for Pioneer 11. The error is de¬ 
termined by use of a five-day batch sequential filter with 
radial acceleration as a stochastic parameter subject to 
white Gaussian noise (~ 500 independent five-day sam¬ 
ples of radial acceleration) [Q||]. No magnitude variation 
of ap with distance was found, within a sensitivity of 2 
x 10 -8 cm/s 2 over a range of 40 to 60 AU. 

Continuing our search for an explanation, we consid¬ 
ered the possibilities i) that the Pioneer 10/11 spacecraft 
had internal systematic properties, undiscovered because 
they are of identical design, and ii) that the acceleration 
was due to some not-understood viscous drag force (pro¬ 
portional to the approximately constant velocity of the 
Pioneers). Both these possibilities could be investigated 
by studying spin-stabilized craft whose spin axes are not 
directed towards the Sun, and whose orbital velocity vec¬ 
tors are far from being radially directed. 

Two candidates were Galileo in its Earth-Jupiter mis¬ 
sion phase and Ulysses in Jupiter-perihelion cruise out of 
the plane of the ecliptic. As well as Doppler, these space¬ 
craft also yielded a considerable quantity of range data. 
Ranging data are generated by cross correlating a phase 
modulated signal with a ground duplicate and noting the 
time delay. Thus, the ranging data are independent of 
the Doppler data, which represent a frequency shift of 
the radio carrier wave without modulation. (For exam¬ 
ple, solar plasma introduces a group delay in the ranging 
data but a phase advance in the Doppler data.) Ranging 
data can be used to distinguish an actual range change 
from a fictitious one caused by a frequency error. 

A quick look at Galileo showed it was impossible to 
separate the solar radiation effect from the anomalous 
constant acceleration with the limited data analyzed (241 
days from 8 January 1994 to 6 September 1994) ||. 

However, an analysis of the radiation pressure on 
Ulysses in its out-of-the-ecliptic journey, from 5.4 AU 
near Jupiter in February 1992 to the perihelion at 1.3 AU 
in February 1995, found a varying profile with distance 
]/]. The orbit solution requires a periodic updating of the 
solar radiation pressure. The radio Doppler and ranging 
data can be fit to the noise level with a time-varying 
solar constant in the fitting model js|]. The inferred so¬ 
lar constant is about 40 percent larger at perihelion (1.3 
AU) than at Jupiter (5.2 AU), a physical impossibility. 
By interpreting this time variation as a true r -2 solar 
pressure plus a constant radial acceleration, we conclude 
that Ulysses was subjected to an unmodelled acceleration 
towards the Sun of (12 ± 3) xlO -8 cm/s 2 . 

With no explanation of this data in hand, our atten¬ 
tion focused on the possibility that there was some er¬ 
ror in JPL’s ODP. To investigate this, an independent 
analysis of the raw data using The Aerospace Corpo¬ 
ration’s Compact High Accuracy Satellite Motion Pro¬ 


gram (CHASMP), which was developed independently 
of JPL’s ODP, was performed. Although by necessity, 
both programs use the same physical principles, plane¬ 
tary ephemeris, and timing and polar motion inputs, the 
algorithms are otherwise quite different. If there were an 
error in either program, they would not agree. (Common 
program elements continue to be investigated.) 

The CHASMP analysis of Pioneer 10 data also showed 
an unmodelled acceleration in a direction along the radial 
toward the Sun ||]. The value is (8.65 ± 0.03) x 10 -8 
cm/s 2 , agreeing with JPL’s result. The smaller error 
here is because the CHASMP analysis used a batch least- 
squares fit over the whole orbit [|J, not looking for a 
variation of the magnitude of ap with distance. 

Without using the apparent acceleration, CHASMP 
shows a steady frequency drift of about —6 x 10 -9 Hz/s, 
or 1.5 Hz over 8 years (one-way only). This equates to a 
clock acceleration, — a t , of — 2.8xl0” ls s/s 2 . Theidentity 
with ap is ap = ate. The drift in the Doppler residuals 
(observed minus computed data) is seen in Figure 1. It 
is clear, definite, and cannot be removed without either 
the added acceleration, ap, or the inclusion in the data 
itself of a frequency drift, i.e., a “clock acceleration” a*. 

If there were a systematic drift in the atomic clocks of 
the DSN or in the time-reference standard signals, this 
would appear like a non-uniformity of time; i.e., all clocks 
would be changing with a constant acceleration. We have 
not yet been able to rule out this possibility. Elements 
common to the Doppler and range tracking systems (e.g., 
DSN station clocks) need to be investigated. For exam¬ 
ple, how and to what accuracy are the clocks at different 
DSN stations tied to each other and to external national 
standards? Are there differences in the orbital fits when 
different stations’ data are analyzed separately? 

Aerospace’s analysis of Galileo data covered the same 
arc as JPL and a second arc from 2 December 1992 to 
24 March 1993. Doppler data from the first arc resulted 
in a determination for ap of ~ (8 ± 3) x 10~ 8 cm/s 2 , a 
value similar to that from Pioneer 10. But the correlation 
with solar pressure was so high (.99) that it is impossible 
to decide whether solar pressure is a contributing factor. 
[Galileo is less sensitive to both the ap- and a t -model 
effects than the Pioneers. Pioneers have a smaller solar 
pressure and a longer light travel time. Sensitivity to a 
clock acceleration is proportional to the light travel time 
squared.] The second arc was 113 days long, starting six 
days prior to the second Earth encounter. This solution 
was also too highly correlated with solar pressure, and 
the data analysis was complicated by many mid-course 
maneuvers. The maneuver uncertainties were so great, a 
standard null result could not be ruled out. 

However, there was an additional result from this sec¬ 
ond arc. This arc was chosen for study because there 
was ranging data. The two-way range change and time 
integrated Doppler are consistent to ~ 4 nr over a time 
interval of one day. This is strong (but not conclusive) 
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evidence that the apparent acceleration is not the result 
of hardware problems at the tracking stations. 

With these added discoveries, what other possible ori¬ 
gins for the signal come to mind? 

One can speculate that there is some unknown interac¬ 
tion of the radio signals with the solar wind. An exper¬ 
imental answer could be given with two different trans¬ 
mission frequencies. Although the main communication 
link on the Ulysses mission is S-up/X-down mode, a small 
fraction of the data is S-up/S-down. We plan to utilize 
this option in further analysis. 

If no normal explanation for the residuals is found, fur¬ 
ther tests of the effect are needed. The weakening Pioneer 
10 signal can still be reacquired for a short time. (The 
NASA/Ames Lunar Prospector Team has intermittently 
done this for training purposes, producing high quality 
data.) Further Ulysses data would also help. 

The Pluto Express mission could provide an excellent 
opportunity for high-quality data from very deep space, 
especially if optical tracking is used. A similar opportu¬ 
nity may exist, out of the plane of the ecliptic, from the 
proposed Solar Probe mission. Under consideration is a 
low-mass module to be ejected during solar fly by. 

With all the above, it is interesting to speculate on the 
possibility that the origin of the anomalous signal is new 
physics 1101. This is true even though the probability is 
that some “standard physics” or some as-yet-unknown 
systematic will be found to explain this “acceleration.” 
This probability is of interest in itself, given that we have 
found no plausible explanation so far. 

The paradigm is obvious. “Is it dark matter or a mod¬ 
ification of gravity?” Unfortunately, neither easily works. 

If the cause is dark matter, it is hard to understand. 
The spherically-symmetric distribution of matter, p ~ 
r -1 , produces a constant acceleration inside the distri¬ 
bution. For this to cause ap , even only out to 50 AU, 
would require the total dark matter to be > 3 x 10” 4 Mq. 
But this is in conflict with the accuracy of the ephemeris, 
which allows only of order a few times 10~ 6 Mq of dark 
matter even within the orbit of Uranus [jllj] . (A 3-cloud 
neutrino model also did not solve the problem |l2[.) 

Contrariwise, the most commonly studied possible 
modification of gravity (at various scales) is an added 
Yukawa force [[b| . Then the gravitational potential is 

V{r) = — GMm[(l + a)r] _1 [l + aexp _r, / A ], (1) 


where a is the new coupling strength relative to Newto¬ 
nian gravity, and A is the new force’s range. Since the 
radial force is F r = —d r V ( r) = ma, the power series for 
the acceleration yields an inverse-square term, no inverse- 
r term, then a constant term. Identifying this last term 
as the Pioneer acceleration yields 

ap = —aia[2(l + a)] -1 |r 2 /A 2 ], (2) 


where a± is the Newtonian acceleration at distance r\ = 1 
AU. (Out to 65 AU there is no observational evidence of 


an r term in the acceleration.) Eq. (||) is the solution 
curve; for example, a = — 1 x 10~ 3 for A = 200 AU. 

It is also of interest to consider Milgrom’s proposed 
modification of gravity 0’ where a oc 1/r 2 for some 
constant ao <C a and a oc 1/r for ao ^ a. Depending on 
the value of the Hubble constant, we find that ao ~ ap. 

Of course, there are (fundamental and deep) theoreti¬ 
cal problems if one has a new force of the phenomenolog¬ 
ical types of those above. Even so, the deep-space data 
piques our curiosity. However, these and other universal- 
gravitational explanations for the Pioneer effect come up 
against a hard experimental wall. 

The anomalous acceleration is too large to have gone 
undetected in planetary orbits, particularly for Earth and 
Mars. NASA’s Viking mission provided radio-ranging 
measurements to an accuracy of about 12 m JT^,|l6[. If 
a planet experiences a small, anomalous, radial accelera¬ 
tion, a a, its orbital radius r is perturbed by 

Ar = -l 6 a A /{GM Q ) 4 -> -r[a A /a N ], (3) 

where l is the orbital angular momentum per unit mass 
and aw is the Newtonian acceleration at r. (The right 
value in Eq. (j^) holds in the circular orbit limit.) 

For Earth and Mars, Ar is about -21 km and -76 km. 
However, the Viking data determines the difference be¬ 
tween the Mars and Earth orbital radii to about a 100 m 
accuracy, and their sum to an accuracy of about 150 m. 
The Pioneer effect is not seen. 

Further, a perturbation in r produces a perturbation 
to the orbital angular frequency of 

Au> = 21 o. a /(GMq) 2&[a A /a,N]- (4) 

The determination of the synodic angular frequency 
uje — is accurate to 7 parts in 10 1 , or to about 
5 ms accuracy in synodic period. The only parameter 
that could possibly mask the spacecraft-determined ap 
is ( GMq ). But a large error here would cause inconsis¬ 
tencies with the overall planetary ephemeris |[TT] . |T7|| . 

We conclude that the Viking ranging data limit any 
unmodelled radial acceleration acting on Earth and Mars 
to no more than 0.1 x 10 -8 cm/s 2 . Consequently, if the 
anomalous radial acceleration acting on spinning space¬ 
craft is gravitational in origin, it is not universal. That 
is, it must affect bodies in the 1000 kg range more than 
bodies of planetary size by a factor of 100 or more. This 
would be a strange violation of the Principle of Equiva¬ 
lence (PE) ju| . The fact an anomalous signal is not seen 
in the analysis of the Viking Lander ranging data gives 
us added confidence that the anomaly is not related to 
DSN hardware. However, the Viking Lander data have 
not been analyzed by either ODP or CHASMP, so we 
cannot make a similar claim regarding software errors. 

Similarly, the A w results rule out the universality of 
the at time-acceleration model. In the age of the uni¬ 
verse, T, one would have ajT 2 /2 ~ 0.7 T. (Another 
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unusual possibility is that there is some unknown non- 
kinematic effect causing a Doppler anomaly.) 

Clearly, more analysis, observation, and theoretical 
work are called for. Further details will appear elsewhere. 
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FIG. 1. Two-way Doppler residuals (observed Doppler ve¬ 
locity minus model Doppler velocity) for Pioneer 10 in mm/s 
vs. time. Solar system gravity is represented by the Sun and 
the planetary systems [JUI. [If one adds one more parame¬ 
ter to the model (a constant radial acceleration) the residuals 
are distributed about zero Doppler velocity with a systematic 
variation ~ 3.0 mm s^ 1 on a time scale ~ 3 months.] The 
outliers on the plot were rejected from the fit. 
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